How does the Eurodollar Interest Rate behave? 
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An empirical analysis on Eurodollar interest rates daily data in the time period 1990-1996, 
is performed and compared with Libor data in the time period 1984-1998. The complementary 
cumulative distributions for the daily fluctuations at different maturity dates and the Power Spectral 
Density are computed. We find that the probability distribution shows 'fat' tails with non-Gaussian 
behaviours. Moreover, we study the correlations among Eurodollar interest rates fluctuations with 
different maturity dates. By using an original clustering linkage, we show how the collective motion 
of the interest rates curve can be analyzed in sub-groups of maturity dates with similar behaviours. 



I. INTRODUCTION 

The last twenty years have seen a wide development of 
studies on the interest rates curves. So far, most of the 
efforts in Finance have been devoted to the study of stock 
prices, where one is dealing with a single object affected 
by stochastic fluctuations. On contrast, the interest rates 
case deals with several objects that follow similar trends 
and can be thought as a curve whose shape is varying in 
time. This makes the subject very challenging since one 
is no more dealing with the statistics of single objects but 
with the motion of a whole complex set. 

Nowadays, these studies are becoming very attrac- 
tive and are approached from many different perspec- 
tives pm]- Several theoretical models have been pro- 
posed [pHyl . These models often assume the fluctuations 
to be Gaussian, neglecting therefore the 'fat tail' effects 
[ 1 2 1 . This is in contrast with the empirical observations 
[ 1 3 1 . The inadequacies of the Gaussian models in finan- 
cial analysis have been reported for a long time p^ , p^[ , 
but the availability of enormous sets of financial data and 
the capability to treat them, has recently renewed the in- 
terest in the subject [p^|-p4[. While the fat-tail property 
of the empirical distribution of stocks price changes has 
been widely documented, the interest rates have been 
less investigated Moreover, the studies of interest 

rates are often limited to only few maturity dates (3 and 
6 months cash rates in J26|, Bund futures in ||). Th° 
paper p7f , is an exception as the US forward rate curve 
with maturity dates up to four years is modelled. 

In the present paper, we perform an empirical analysis 
of time series based on daily prices of Eurodollar future 
contracts in the time period 1990-1996. The data are 
described in Section [n| (more details can be found in the 
paper of Bouchaud et al. |2^]). Our study is mainly fo- 
cused on the 'tail' regions of the fluctuations probability 
distributions and on the correlations among time series 



at different maturity dates. The results concerning the 
probability distribution and the Power Spectral Density 
of interest rates fluctuations are presented in Section [II 
and in Section IV, respectively. In Section the cor- 
relations are analyzed in terms of an original clustering 
linkage. Finally, some remarks are reported in the last 
Section. 



II. EMPIRICAL DATA 

An Eurodollar is a dollar deposited in a bank out- 
side the United States by another bank. The Eurodol- 
lar future contract is a future contract on an interest 
rate. These contracts are traded on an exchange. The 
three-month Eurodollar futures traded on the Chicago 
Merchantile Exchange (CME) is the most popular of 
the futures contract on short-term interest rates. These 
contracts trade with delivery months of March, June, 
September, and December up to ten years into the future. 
The variable underlying the contract is the Eurodollar in- 
terest rate applicable to a 90-days period beginning on 
the third Wednesday of the delivery month. When the 
delivery date is reached, and the actual interest rate for 
the 90-day period is known, the contract is settled in 
cash. The three-month Eurodollar interest rate is the 
same as the three-month Libor (London Interbank Of- 
fered Rate), that is the rate at which large international 
banks found much of their activities. Specifically, Libor 
is the rate at which one large bank is willing to lend 
money to another large bank |l|,^,^8|,^9| . All these rates 
are determined in trading between banks and change as 
economic conditions change. 

Following the line of Bouchaud et al. p7| , [30|]3"l| ] , we con- 
sider three-months future rates as instantaneous forward 
rates so obtaining daily values for the forward interest 
rate curve f(0, t), where t is the current date and 9 is the 
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maturity date. In the time period analyzed 1990-1996, 
the variable t counts incrementally the trading days be- 
tween 1 and T = 1831. On this time period we have 
16 different time series corresponding to several matu- 
rity dates ranging from 9 = 3 to 48 months with a step 
of three months. The behaviour of f{9, t) as function of t 
is shown in FigJ^ where, for a better visualisation of the 
plot, we report only those values corresponding to the 
following maturity values: 9 — 3, 15,30,48 months. The 
interest rates behaviours for the other maturity dates fol- 
low very similar trends in time, and stay mostly inside 
the shape traced by the two extreme maturity values, 
namely 9 = 3 and 9 = 48 months. 



III. PROBABILITY DISTRIBUTIONS 

To investigate the nature of the stochastic process un- 
derlying the interest rates changes, we first analyze the 
probability distribution of the daily changes in the inter- 
est rates. 

Let us start by the computation of the standard devi- 
ation of Af(9,t), which is defined as: 
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FIG. 1. Eurodollar interest rates (f(0,t)) as function of t 
for 9 = 3, 15, 30, 48 months. 

The interest rate fluctuations are analyzed by studying 
the changes from one day to the following day: 



Af(9,t) = f(9,t + At)-f(9,t) 



(1) 



where At = 1 day. In Fig.|2j these fluctuations are re- 
ported as a function of t for f(9 = 3,t). All the others 
at different maturity dates have similar behaviours char- 
acterised by stochastic fluctuations around the zero. 
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where T\ and T2 delimit the range of t, and < A/ > is 
the average over time of Af{9, t) (which tends to zero for 
T 2 - Tx -> 00). In the whole period 1990-1996, a{9) v.s. 
9 has a minimum value of 0.063 at 9 = 3 months, then 
increases up to a maximum value of 0.087 at 9 = 1 year, 
once reached this maximum the value slowly decreases 
until a (6) = 0.062 at 6 = 4 years. 
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FIG. 3. Standard deviation a(9) as function of the years, 
for Eurodollar at 9 = 3, 6, 15, 30, 48 months (time period 
1990-1996) and Libor for 9 = 1, 3, 6 months (time period 
1984-1998). 

We compare the Eurodollar standard deviation of A/ 
for each single year between 1990-1996 with the one of the 
Libor data corresponding to a longer time period from 
1984 to 1998. In Fig.| the values of a{9) in the different 
years are plotted for several maturity dates. The stan- 
dard deviations for both Eurodollar and Libor series with 
the same maturity dates 9 =3 and 6 months, have similar 
values and show decreasing trends between 1984 to 1998. 

The statistical distributions of the daily changes in 
the interest rates are better analyzed by computing their 
complementary cumulative distributions ^(A/), that tell 
us the probability to find a daily change which is larger 
than Af: 



FIG. 2. Eurodollar interest rate fluctuation (Af(9,t)) as 
function of t for 9 = 3 months. 
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with being the probability density distributions of 
Af(0,t). 
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FIG. 4. Complementary cumulative distribution >]/(A/) 
v.s. Af for Eurodollar data at different maturity dates: 9 = 3, 
15, 30, 48 months. 

We find that our time series are highly leptokurtic 
for all the maturity dates and are characterized by non- 
Gaussian profiles for large interest rates changes. The 
results for 'J? (A/) are reported in Fig.^ for the tails re- 
gions for positive changes of Af with 8 = 3, 15, 30, 48 
months. These curves show, in a log-log scale, three dis- 
tinct parts characterised by different decreasing trends: 
i) a first region (0.01 < |A/| < 0.05) with slopes be- 
tween —0.4 -. — 0.6; ii) an intermediate region (0.05 < 
|A/| < 0.2) with a faster decreasing behaviour (slopes 
-1.6 -2.7); iii) a final region (|A/| > 0.2) which is 
characterized by larger slopes (between -2.3 up to -10). 
Comparable results are obtained for all the time series at 
different maturity dates. 



In Fig.|[ we plot the tail region of *(A/) at 9 = 3 
months as function of both positive and negative changes 
of Af. We observe that the probability distributions 
are slightly asymmetric with negative changes that are 
slightly more probable than the positive ones. On the 
other hand, the large variations result more probable for 
the positive changes, and the overall distribution leads 
-correctly- to a zero mean value for Af. In the same fig- 
ure the two empirical curves for the Eurodollars are com- 
pared with the two curves for the Libor at = 3 months. 
Large fluctuations in Eurodollar and Libor data show 
very similar behaviours with comparable probability val- 
ues. In these distributions the tails are non-Gaussian: 
they are fatter than the Gaussian ones. 

Instead that with the Gaussian distribution, a more 
appropriate comparison should be performed with the 
general class of Levy, Khinthine stable distributions |32| . 
A symmetric stable distribution with zero mean can be 
written as [01: 
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where < a < 2 and 7 is a positive scale factor. For 
a = 2 and 7 = er 2 /2, Pj, is a Gaussian distribution with 
zero mean and standard deviation equal to a; whereas 
for a = 1, Pi, is the Lorentzian distribution. 

In Fig.|^ we compare the Eurodollar empirical data 
with a symmetric stable distribution, given by Eq.(|]), 
having a — 1.45. One can note that these curves have 
a tail region that seems to follow a Levy stable distri- 
bution. But we may expect a faster decrease to zero 
for larger fluctuations of Af, as for instance the one de- 
scribed by the truncated Levy distribution ]33| . A deeper 
analysis on this aspect will be done in a future work. Let 
note that all time series at different maturity dates show 
behaviours analogous to the one at 9 — 3 months. 
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FIG. 5. Complementary cumulative distribution (>3/(A/)) 
and cumulative distribution (1-^(A/)) v.s. |A/|, for Eurodol- 
lar and Libor data at 9 = 3 months. The grey-line refers to 
the Gaussian behaviour, whereas the black-line is the Levy 
stable distribution. 
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FIG. 6. Complementary cumulative distribution ^(A/) 
v.s. Af for Libor data at the maturity date 9 — 3 months, 
for different values of the time increment At. 
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Fig.|| reports ^(A/) for the Libor non-overlapping Af 
time series with 8 = 3 month and At values ranging be- 
tween 1 and 10 (see Eq.([j])). Note that the Gaussian 
behaviour tends to appear in the tail regions when the 
At increases. 



IV. SPECTRUM ANALYSIS 



In this section we investigate the statistical properties 
of our time series in the frequency domain. To this end 
we compute the power spectral density (PSD) Q] of our 
random variables by using the periodogram approach to 
PSD estimation, that is currently one of the most pop- 
ular and computationally efficient PSD estimator. This 
periodogram spectral estimate is obtained as the squared 
magnitude of the output values from an FFT performed 
directly on the data set. The power spectrum analysis 
for f{9,t) indicates a stocastic signal with spectral com- 
ponents which are decreasing as a power law: 



(5) 



with f3 = 1.8. Fig.^ shows the power spectrum of the spot 
interest rate r(t) = f{9 = 3,t) on a log-log scale. The 
slope of the linear fit is j3 = —1.77 ± 0.01 and, within 
the error, it is the same for all the other Eurodollar in- 
terest rates at different maturity values, and it is in good 
agreement with the Libor at 9 = 1, 3, and 6 months 
which have = -1.80 ± 0.01. 
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FIG. 7. Power spectrum S(uj) of the Eurodollar interest 
rates at 9 — 3 months, as a function of the frequency ui. 

The same analysis for Af(9,t) (see Figj|) shows a flat 
spectrum, from which we can infer that the interest rates 
changes are "pair wise uncorrelated" and the stochastic 
process is approximately a white noise. 
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FIG. 8. Power spectrum S(u>) of the Eurodollar inter- 
est rates fluctuations at 6 = 3 months, as a function of the 
frequency u). 



V. CORRELATIONS 

In order to investigate the correlations among the in- 
terest rates daily changes at different maturity dates, let 
us here define the cross-correlation matrix: 



< AfiAfj >-<Af l >< Af 3 > 



(7,(7, 



(G) 



where i,j = 1, .... 16, Afi indicates the interest rate daily 
changes for maturity 9i = 3i months, and at is the stan- 
dard deviation of interest rates changes Afi. The sym- 
bol < ... > denotes a time average performed over all the 
days of the investigated period. The correlation coeffi- 
cients are computed between all the pairs of indices la- 
belling our Eurodollar series. Therefore we have a 16 x 16 
symmetric matrix with Cj $ = 1 on the diagonal. 

The daily changes in the interest rates are highly cor- 
related with coefficients in the range 0.48 < Cij < 0.98 
(for i 7= j). These strong correlations between the inter- 
est rates fluctuations suggest that one can consider their 
behaviour as a collective motion of a set of quantities 
joined together with some interaction force. 
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FIG. 9. Correlation coefficients (cij, with i,j = 1..16 and 
i<j) for Af{6,t) and/(M). 
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The correlation coefficients Cij for A/j = Af(9i,t) are 
shown in FigJ| together with the correlation coefficients 
for f(9i,t) (which are obtained substituting A/(0j,t) 
with f(9i,t) in Eq.(||)). The plot reads in the following 
way: the vertical lines group together the correlation co- 
efficients Ci.j for a fixed value of i and j ranging between 
i + 1 to 16. Looking at the plot from left to right, i in- 
creases its value from i = luptoi = 15. The coefficients 
Cjj- have decreasing trends with j. This corresponds to 
the intuitive fact that interest rates with close maturity 
dates are more correlated than interest rates with distant 
maturity (e.g. c% t 2 > C1.4 > ci.% > ci.iq). Such decreas- 
ing trends support the idea that interest rates are di- 
rectly interacting between neighbours only (which means 
that the maturity 9i has direct interactions with maturity 
9i±i), so making a sort of 'string' 
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A closer look at Fig.|] shows that the decreasing trend 
in the correlation coefficients for A/j, is not perfectly 
monotonia In particular, one can observe deviations 
from the monotonic trend for correlation coefficients in- 
volving j > 8 (i.e. maturity dates 9j > 2 years). 

To understand the geometrical and topological struc- 
ture of the correlation coefficients, we introduce the met- 
ric distance dij |3^] between the series A/, and Afj 
which is defined as: 



2(1 



(7) 



By definition, Qj- is equal to zero if the interest rates 
series i and j are totally uncorrelated, whereas Cjj = ±1 
in the case of perfect correlation/anti-correlation. There- 
fore, 



di,j =0 if Cjj = 1 
d it j = V2 if 
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Note that this metric distance fulfils the three axioms of 
a metric: 1) dij — if and only if i = j; 2) dij = dj.f, 
3) dij < di t k + dk,j- Indeed, dij is the Euclidean dis- 
tance, in a T-dimensional space, between the two vectors 

with components [Af(9 t , l),A/(0 i5 2), , Af(9 ll T)]/a i 

and [Af(9j, l),....,A/(6*j, T)]/aj. These vectors have uni- 
tary lengths. Therefore, they connect the centre of a unit 
sphere, in a T-dimensional space, with points on its sur- 
face. Similar behaviours of the interest rate changes cor- 
respond to two close points on the sphere (dij « 1). Un- 
correlated behaviours correspond to two vectors pointing 
toward perpendicular directions (dij = \/2). Whereas 
anti-correlated behaviours are associated to two vectors 
pointing toward two opposite directions (dij = 2= sphere 
diameter). 

We study the geometrical and topological arrange- 
ments of Eurodollar interest rates by introducing the 
ultra-metric distance dij which satisfies the first two 
properties of the metric distance and replaces the tri- 
angular inequality with the stronger condition: dij < 



maxldi^jdk.j], called 'ultra-metric inequality'. Once the 
metric distance dij between our interest rates fluctua- 
tions is known, one can introduce several ultra-metric 
distances. 

Mantegna et. al have used the 'subdominant ultra- 
metric', obtained by calculating the minimum spanning 
tree [|tJ connecting the stock indices [f36|,[38j . In this pa- 
per we consider a different ultra-metric space that em- 
phasizes the cluster-structure of the data. 

In our case, a "cluster" is a set of elements with rela- 
tive distances dij which are smaller than a given thresh- 
old distance 5, whereas disjoined clusters have some ele- 
ments which are at distances larger than 6. Let us define 
the size of a cluster as the number of elements that it 
contains, and its diameter as the maximum relative dis- 
tance between them. Therefore, an isolated element can 
be considered a cluster of size 1 and zero diameter. 

We define the ultra-metric distance dij between two 
elements i,j belonging to two different clusters as the 
maximum metric distance between all the couples of el- 
ements in the two clusters. Therefore, the ultra-metric 
distance between two isolated elements (clusters of size 1) 
coincides with the metric distance. Whereas, the ultra- 
metric distance between two elements inside the same 
cluster is equal to the cluster diameter. 

The linkage procedure that we adopt is described in 
details in Appendix [A|. This procedure yields to a set of 
clusters which have diameters which are smaller or equal 
than the threshold distance 5, and have relative ultra- 
metric distances d which are larger than S. 
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FIG. 10. Hierarchical tree obtained from the correlation co- 
efficients of the 16 Eurodollar interest rates fluctuations time 
series in the time period 1990-1996. (On the a;-axis are re- 
ported the maturity dates and on the j/-axis the ultra-metric 
distances.) 

The result of this recursive linkage procedure, for the 
16 maturity dates (9 = 3, ..,48) in the whole time pe- 
riod 1990-1996, is reported in the hierarchical tree shown 
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in Fig.[10[ The clustering begins with maturities 39-45 
that gather together at 8 = 0.185; at a little larger value 
6 = 0.186, the couple 42-48 makes another cluster and the 
couple 30-36 merges at the 8= 0.206; then the maturity 
33 joins the first cluster at a 8 — 0.216. The clustering 
process goes on (see Fig.|l0|) and ends when only one large 
cluster is formed at 6 — 1.02. 
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FIG. 11. Grey-scale-tones visualisation of the clusters at 
different 8 values, for the 16 Eurodollar interest rates fluctu- 
ations time series in the time period 1990-1996. 

This mechanism of nucleation, growth, and coalescence 
is better visualised in Fig. pd| where the clusters are explic- 
itly drawn with different grey-tones for 6 distinct ranges 
of 5 (0-7- 0.2; 0.2 -j- 0.4; 0.4 -j- 0.6; 0.6 4- 0.8; 0.8 -h 1; 
1 -4 1.2). It is evident from Figs Jl^ and |ll| that, for ultra- 
metric distances 8 in the range between 0.4 and 0.6, the 
data set is gathered into 3 main clusters: Cls\ — {3,6}, 
Cls 2 = {9,12,15,18,21}, Cls 3 = {24, 27, 45, 48}, 
with diameters d = 0.553, 0.411, 0.433, respectively. The 
first cluster {Cls\) gathers together Af(0,t) with matu- 
rity shorter than 1 year; Cls2 contains those with matu- 
rity dates between 1 year and 2 years; whereas Cls^ has 
those with maturity dates which are larger than 2 years. 
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FIG. 12. The nucleation distance 8 n and the distance S p 
as functions of the years between 1990-1996. 



Let us now proceed looking at the clusters arrangement 
of Af(6, t) for each single year between 1990 to 1996. We 
find that the cluster arrangement in each single year is 
similar to the one observed in Figs, 10, ll] where the whole 
time period is considered. First we have the nucleation of 
a couple of maturity dates at distances 8 = 8 n whose val- 
ues are between 0.109 and 0.210 (see Fig.[l2], where these 
nucleation distances are plotted for each year). Then we 
observe the growth of several clusters, that in the dis- 
tances range between 0.3 to 0.6, makes a structure of 
four well identifiable clusters. The first one is composed 
by the maturity 6 — 3 and sometime 9 — 6. The second 
includes time series with maturity date values around 1 
year (9 — 6 4 15). The third gathers maturity dates 
between 2 and 2.5 years (0 — 12 33) and the fourth 
the largest maturity dates > 2 years (8 = 27 -j- 48). Fi- 
nally, these clusters coalesce into a single one at distances 
8 = S p , whose values (ranging between 0.852 and 1.131) 
are plotted in Fig jT^ for each year. Note that the val- 
ues of 8 n and 8 P , calculated for the whole time period 
(1990-1996), are within the range of the values plotted in 
FigM 



VI. CONCLUSIONS 



We have empirically analyzed Eurodollar interest rates 
at different maturity dates in the time period 1990-1996. 
We observed that the tail region has an evident non- 
Gaussian behaviour - a crucial feature for the risk man- 
agement analysis. Therefore, by using a Gaussian prob- 
ability density function one underestimates the proba- 
bility of large interest rates fluctuations. These results 
are confirmed by the same analysis done on Libor time 
series in the longer time period from 1984 to 1999. More- 
over, a power spectrum analysis has been carried out in 
order to give some insights into the stochastic process 
underlying the examined interest rates time series. The 
power spectrum for the interest rates fluctuations shows a 
white-noise-like flat behaviour. Whereas, the one for the 
interest rates has a power law spectrum with exponent 
1.77 (Eurodollar), and 1.80 (Libor) for all the different 
maturity dates. 

We have also analyzed the correlations between inter- 
est rates fluctuations at different maturity dates by using 
a linkage procedure that emphasizes the cluster struc- 
ture. In the time period 1990-1996, we have found the 
formation of three clusters joining together time series 
corresponding to three maturity dates periods: shorter 
than 1 year, between 1 to 2 years, larger than 2 years. A 
cluster structure with similar features have been found 
for each single year. The high correlation values that we 
have found, confirm that the interest rates fluctuations 
behaviour at different maturity dates must be seen as a 
collective motion. On the other hand, the result from 
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the clustering linkage procedure indicates that the de- 
viations from this collective behaviour can be gathered 
in some maturity-groups. This might suggest that the- 
oretical models for the interest rates behaviours should 
usefully consider different kinds of interactions in differ- 
ent sub-groups of maturity dates. 
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APPENDIX A: LINKAGE PROCEDURE 



At the starting point of our linkage procedure we have 
a set of N(= 16) isolated elements which can be consid- 
ered as N clusters of size equal to 1 and diameter equal 
to zero. We define the ultra-metric distance dij between 
two elements i,j belonging to two different clusters as 
the maximum metric distance between all the couples of 
elements in two clusters. Whereas, the ultra-metric dis- 
tance between two elements inside the same cluster is 
equal to the cluster diameter. 

The first step of the linkage procedure finds the couple 
of elements i,j which are at the smallest distance dij. 
This couple is linked together in a new cluster of size 
2 and diameter dij. Then we seek for the pair of ele- 
ments with the next-smallest ultra-metric distance, and 
we generate a new cluster, or we enlarge the previous one 
if one of the two clusters belongs to that cluster. This 
procedure is recursively iterated for increasingly larger 
distances by joining isolated elements to clusters or by 
merging together larger clusters. The procedure termi- 
nates when all the elements are jointed together in a sole 
cluster. This ending diameter coincides with the largest 
metric distance between two elements in the data set. 
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